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ABSTRACT 

The recent detection of significant linear polarization at mm and sub-mm wavelengths in 
the spectrum of Sgr A* (if confirmed) will be a useful probe of the conditions within several 
Schwarzschild radii (r\s) of the event horizon at the Galactic Center. Hydrodynamic simulations 
of gas flowing in the vicinity of this object suggest that the infalling gas circularizes when it 
approaches within 5 — 25 r$ of the black hole. We suggest that the sub-mm "excess" of emission 
seen in the spectrum of Sgr A* may be associated with radiation produced within the inner 
Keplerian region and that the observed polarization characteristics provide direct evidence for 
this phenomenon. The overall spectrum from this region, including the high-energy component 
due to bremsstrahlung and inverse Compton scattering processes, is at or below the recent 
Chandra measurement, and may account for the X-ray source if it turns out to be the actual 
counterpart to Sgr A*. 

Subject headings: accretion — black hole physics — hydrodynamics — Galaxy: center — magnetic 
fields: dynamo — radiation: polarization 
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1. INTRODUCTION 

Discovered over 25 years ago (Balick & Brown 1974), Sgr A* is a bright, compact radio source coincident 
with the dynamical center of the Galaxy, and now provides possibly the most compelling evidence for the 
existence of supermassive black holes. The suggested central dark mass concentration within the inner 
0.015 pc of the Galactic center is 2.6 ± 0.2 x 10 6 M (Genzel et al. 1996; Eckart & Gcnzcl 1996; Eckart & 
Genzel 1997; Ghez et al. 1998). (0.1" corresponds to 800 Astronomical Units, or roughly 1.2 x 10 16 cm at a 
distance of 8.5 kpc.) Most of this mass is probably associated with Sgr A*. 

The spectrum of this unusual object is seen to be bumpy, but can be described as a power-law with 
a spectral index a that varies between roughly 0.19 — 0.34 (S u cx v a ) between GHz and mm wavelengths. 
However, one of the most interesting features currently under focus is the suggestion of a sub-millimeter 
(sub-mm) bump in the spectrum (Zylka et al. 1992; Zylka et al. 1995), since in all emission models the 
highest frequencies correspond to the smallest spatial scales, so that the sub-millimeter emission comes 
directly from the vicinity of the black hole (Mclia 1992, 1994; Mclia, Jokipii & Narayanan 1992; Coker and 
Mclia 2000). The existence of this bump (or "excess") has been uncertain due to the variability of Sgr A*, 
but is now well established following a set of simultaneous observations (from A20cm to Almm) using the 
VLA, BIMA, Nobcyama 45 m, & IRAM 30 m telescopes (Falcke, et al. 1998). 

More recently, radio observations of Sgr A* have focused on the detection of polarization from this 
source. Although the upper limits to the linear polarization in Sgr A* are found to be quite low (less than 
1%) below 86 GHz (Bower et al. 1999), this is not the case at 750, 850, 1350, and 2000 ^m, where a 
surprisingly large intrinsic polarization of over 10% has now been reported (Aitken, et al. 2000). ^From 
the lack of polarization at longer wavelengths, Aitken et al. conclude that their measured values at higher 
frequencies must arise in the mm/sub- mm "excess". These observations also point to the tantalizing result 
that the position angle changes considerably (by about 80°) between the mm and the sub-mm portions of 
the spectrum, which one would think must surely have something to do with the fact that the emitting gas 
becomes transparent at sub-mm wavelengths (Melia 1992, 1994). 

In a companion paper (Melia, Liu & Coker 2001) we suggested that the mm and sub-mm "excess" 
in the spectrum of Sgr A* may be the first indirect evidence for the anticipated circularization of the gas 
falling into the black hole at 5 — 25 r$, where r$ = 2GM/c 2 is the Schwarzschild radius. The abundance 
of gas in the environment surrounding Sgr A* clearly points to accretion as the incipient cause of the 
ensuing energetic behavior of this source (Melia 1994) — whether or not it eventually leads to expulsion of 
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some plasma at smaller radii (see, e.g., Falcke et al. 1993). In their simulation of the Bondi-Hoyle accretion 
onto Sgr A* from the surrounding winds, Coker & Melia (1997) concluded the accreted specific angular 
momentum I = Xrsc can vary by 50% over 200 years with an average equilibrium value in A of about 
30 or less. Although this shows that relatively little specific angular momentum is accreted — so that large 
disks (such as those required in ADAF models; Narayan et al. 1996) probably do not form around Sgr 
A* — it does nonetheless lead to the expectation that the plasma must circularize toward smaller radii 
before flowing through the event horizon. However, given the fluctuations in the accreted value of A (both 
in magnitude and sign!), this Kcplcrian flow is variable, and it probably dissolves and reforms (possibly 
with a different sense of spin) on a time scale of <~ 100 years or less. 

Melia et al. (2001) demonstrated how this dichotomy comprising a quasi-spherical flow at radii beyond 
50 rs or so, and a Keplerian structure toward smaller radii, may be the explanation for Sgr A*'s spectrum, 
including the appearance of the "excess", which is viewed as arising primarily within the circularized 
component. It is our intention in this Letter to demonstrate how the linear polarization data may now be 
taken as direct evidence for the existence of this accretion profile. 



The structure of the flow within the circularization radius (at ~ 5 — 50 rs) is developed fully in 
Melia et al. (2001). For the sake of completeness, we highlight the key elements of this calculation 
below. Central to the modeling of the sub-mm "excess" is the supposition that within the Keplerian flow, 
a magnetohydrodynamic dynamo produces an enhanced (though still sub-equipartition) magnetic field, 
dominated by its azimuthal component (Hawley, Gammie & Balbus 1996). Although the process of magnetic 
field dissipation suppresses the field intensity well below its equipartition value in the quasi-spherical region 
(Kowalenko & Melia 2000), the magnetic dynamo evidently overwhelms the rate of field destruction in 
the differentially rotating portion of the inflow, and the field reaches a saturated intensity. The sub-mm 
"excess" in Sgr A* may be the thermal synchrotron radiation produced in this inner region. 

At a (cylindrical) radius r in the Keplerian flow where the column density is S and the angular velocity 
is fl = (GM/r 3 ) 1 / 2 , the radial velocity is given as (e.g., Stoeger 1980) 



where v = (2/3)W ri p/Y> Q is the kinematic viscosity, and W ri p is the vertically integrated sum of the Maxwell 
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METHODOLOGY 




(1) 
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and Reynolds stresses (Balbus et al. 1994). For the problem at hand, the Maxwell stress dominates, and 

Wr+ « A, J dz (^-) , (2) 

where B is the turbulent magnetic field (the average inside the integral being taken over time). Numerical 
simulations (e.g., by Brandenburg et al. 1995) show that f3 v changes very slowly with r. In the particular 
cases considered by these authors, (3 V ranged in value from w 0.1 to 0.2, while r decreased by a factor of 5. 
For simplicity, we will here adopt a "mean" value of ~ 0.15 for this quantity. 

For steady conditions, one can obtain the vertical profile by assuming that the gas is in local hydrostatic 
equilibrium. Balancing gravity and the pressure gradient in the vertical direction, we obtain the scale height 
H = ^/2R g Tr 3 / fiGM, where T is the gas temperature at radius r, R g is the gas constant, and /x is the 
molecular weight. For simplicity, we will assume that the Kcplcrian flow is axisymmetric and is independent 
of the vertical coordinate. Written another way, we have (HVl) 2 = 2P/p, where P is the gas pressure and 
p is the mass density of the gas. The numerical simulations of the magnetohydrodynamic dynamo effect 
indicate that the field intensity is somewhat below its equipartition value, so that 

J dz{—)*(3 v J Pdz = (3 p ^L, (3) 

where [3 P is roughly constant with a value of w 0.02. Thus, with M = — 2-7rr £iv, we can integrate Equation 
(111) to obtain v r , and T follows directly from the energy conservation equation (Melia et al 2001). 



3. CALCULATION OF THE SPECTRUM 

The flux density (at earth) produced by the Keplerian portion of the flow is given by 



F V0 = JpJ I v V^-rs/r dA , (4) 



where D = 8.5 kpc is the distance to the Galactic Center, v$ is the observed frequency at infinity and v' 
is the frequency measured by a stationary observer in the Schwarzschild frame. (For simplicity, we here 
assume the metric for a non-spinning black hole. A more thorough exploration of the parameter values, 
including the black hole spin, will be discussed elsewhere.) The frequency transformations are given by 



1-vl/c* 



vo = v'^l - r s /r , v' = v- — -, (5) 

1 - [Vj /cj cos 

where v is the frequency measured in the co-moving frame, and 9 is the angle between the velocity 
and the line of sight. Since the radial velocity is always much smaller than v$, we ignore this component 
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in the transformation equations. So cos 6 = sini cos^>, where i is the inclination angle of the axis 
perpendicular to the Keplerian flow, and <j> is the azimuth of the emitting element. When the Doppler 
shift is included, the blue shifted region is located primarily near = while the red shifted region is at 
4> = 7r. The other quantities that are necessary for an evaluation of the flux density are the area element 
dA = yjl — rs/r 1 cosi r dr d(f>, and the specific intensity I v > = B' V ,{1 — e~ T ), where 



K> = i \.. ,l^a B v , (6) 



and the optical depth is 



1 — (^0 /c) cos 9 



ds = Kv ^ 1 -^ Cme , (7) 



where k v is the absorption coefficient. When r <C 1, Kirchoff's law allows us to write 



I v , « B' v ,r = e v —. I v ^ | , (8) 

cos i y 1 — (v^/c) COS ' 

where e v = B u k v is the emissivity. 

The presence of a substantial azimuthal component of the magnetic field makes it convenient to 
calculate the observed flux directly from the Extraordinary and Ordinary components of the intensity. The 
most convenient approach is to select the symmetry axis of the Keplerian flow as the reference direction. 
The observed flux densities in the azimuthal and the reference directions are given by 

F lva = Jpf (# cos 2 (/>' + I° v , sin 2 4>')y/l - rs/r dA , (9) 

F 2vo = ^2 / sin2 <f>' + C cos 2 0V1 - rs/r dA , (10) 

respectively, where </>' + tt/2 is the position angle of the magnetic field vector within the emitting element 
that has an azimuth of </>, so that cot0' = cot</> cosi. If,, and 1°, are the specific intensities for the 
Extraordinary and Ordinary waves, respectively. For thermal synchrotron radiation, the emissivities are 

V3e 3 



e = 



7 Bsm6' N(E)[F(x) + G(x)] dE , (11) 
Jo 



&Ttm e c 2 J 

e° = V fisinfl' / N(E)[F(x) - G(x)} dE , (12) 
8Trm e c z J 

where N(E) is the electron distribution function at energy E, and 

, cosO — v^/c Airvmlc 5 . . 

C ° = l-(v^/c)cos9 ' X= 3eBsm9>E 2 ' ( U > 



pOO 

F(x) =x K 5/3 (z) dz , G(x) = x K 2/3 (x) 

J X 



(14) 
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fT 5 /3 an d ^2/3 are the corresponding modified Bessel functions (Pacholczyk 1970). The total flux density 
produced by the Keplerian portion of the flow is the sum of these two. The expected fractional polarization 
is then given by P Vo = (F Xva - F 2vo )/(F ivo + F 2uo ). 

The temperature in this region reaches ~ 10 11 K, for which inverse Compton processes must be taken 
into account. The self-Comptonization of the sub-mm radiation is calculated according to the prescription 
in Melia et al. (2001), based on the algorithm described in Melia & Fatuzzo (1989). 



4. RESULTS AND CONCLUSIONS 

The best-fit model for the polarized mm and sub-mm emission from Sgr A* (Aitken et al. 2000) is 
shown in Figures 1 (the inset) and (the solid curve of) 2. The peak frequency of the flux density is 2.4 x 10 11 
Hz, and the flip frequency (at which the position angle changes by 90°) is 2.8 x 10 11 Hz. Below this 
frequency, the first component is smaller than the second, and the corresponding percentage polarization is 
therefore (by definition) negative. Above the flip frequency, the first component is larger. Although the fit 
is not optimized, both the spectrum and the percentage polarization appear to be consistent with the data. 
It is to be noted that the peak frequency is actually smaller than the flip frequency, which is distinct from 
other models that may also produce a rotation of the position angle (see Aitken et al. 2000). 

It is rather straightforward to understand the polarization characteristics in this model. In the optically 
thick region (below about 1.6 x 10 11 Hz), the specific intensity of the Extraordinary and Ordinary waves is 
almost isotropic in the co-moving frame because the optical depth r is very large. Even with the inclusion 
of the Doppler effect, the emissivity of the source is relatively independent of position angle. But the 
optical depths are different for the two waves, as indicated by Equation (Q), and the specific intensity of the 



Extraordinary wave is slightly larger than that of the Ordinary wave. From Equations (11) and (|12|), we 
see that the second component is larger than the first, and the percentage polarization is therefore negative 
according to the definition of P„ . With an increase in frequency, the Extraordinary amplitude becomes 
even larger (relative to that of the Ordinary wave) and so the percentage polarization increases. 

However, in the optically thin region, the specific intensity is given by Equation (|J). The synchrotron 
emissivity is very sensitive to the angle between the line of sight and the magnetic field vector B; synchrotron 
radiation is beamed into a plane perpendicular to B in the co-moving frame. With the inclusion of the 
Doppler effect, the radiation is beamed into a cone, and the dominant contribution comes from the blue 
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shifted region which has an azimuth of about zero. Therefore, since the Extraordinary wave is more intense 
than the Ordinary wave and the integrals ( |Tl| ) and ( |l2| ) are dominated by radiation from the emitting 
element with an azimuth of about zero, the first component is larger than the second. In this case, the 
fractional polarization becomes positive. 

In other words, the optically thick emission is dominated by emitting elements on the near and far sides 
of the black hole, for which the Extraordinary wave has a polarization direction parallel to the reference 
axis. In contrast, the dominant contribution in the thin region comes from the blue shifted emitter to the 
side of the black hole, where the Extraordinary wave has a polarization direction mostly perpendicular 
to this axis. The sharp decrease in polarization at still higher frequencies is due to the diluting effects 
of Comptonization emission which begins to dominate over Synchrotron emission at that point. The 
inclination angle dependence of the fractional polarization associated with emission by the Keplerian 
portion of the inflow is shown in Figure 2. 

Several issues remain to be investigated. An important result of our analysis is that only modest 
accretion rates appear to be consistent with the polarization characteristics of Sgr A* at mm and sub-mm 
wavelengths. The emitting region is compact — evidently no larger than a handful of Schwarzschild radii. 
Yet hydrodynamical simulations (Coker & Melia 1997) suggest that the rate at which plasma is captured 
at larger radii (of order 10 r$ or so) is several orders of magnitude higher. If our modeling is correct, this 
would seem to suggest that M is variable, perhaps due to a gradual loss of mass with decreasing radius 
(see, e.g., Blandford & Begelman 1999). It is essential to self-consistently match the conditions within the 
Keplerian region of the flow with the quasi-spherical infall further out. These calculations are currently 
under way, and the results will be reported elsewhere. In addition, if the Chandra source is indeed the 
counterpart to Sgr A* (Baganoff et al 2000), then the spectrum shown in Figure 1 suggests a correlated 
variability between the sub-mm and X-ray fluxes, which can be tested with the next round of coordinated 
observations. 
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Fig. 1. — The overall spectrum from this region (optimized to fit the Nobeyama and IRAM data; Falcke et al. 
1998) is shown in the main portion of this figure. Dashed: thermal synchrotron, dotted: self-Comptonization, 
solid: total spectrum. The bremsstrahlung emission is negligible on this scale. The parameter values are 
M = 1 x 10 16 g s" 1 , P p = 0.036, (3 V = 0.27, r t = 1.0 r s and r = 5.0 r s . The inclination angle of the axis 
perpendicular to the Keplerian plane is i = 60°. ft is also necessary to specify the ratio of v r to its free-fall 
value at tq. For this model, this ratio is 5.0 x 10~ 5 . The mm to sub-mm spectrum corresponding to the best 
fit model for the JCMT data (which were obtained 3 years after the Nobeyama plus IRAM observations) is 
shown in the inset. The dotted curve corresponds to the first component and the dashed curve corresponds 
to the second component. The solid curve is the sum of these two. Here, M = 4.0 x 10 15 g s _1 , (3 P — 0.02, 
/?„ = 0.15, n — 1.0 rs, r = 4.0 r s and i = 45°. The ratio of velocities at r is 4.0 x 10~ 5 . 
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Fig. 2. — The percentage polarization for the best fit model whose spectrum is shown in the inset of Fig. 1 is 
here indicated by the solid curve (i — 45°). The inclination angle dependence of this quantity is illustrated 
by the curves corresponding to an inclination angle i = 35° (dotted) , and i = 55° (dashed) . 



